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Abstract 

The fermionic Ising spin glass models in a transverse field are inves- 
tigated in a Grassmann path integral formalism. The Parisi's scheme of 
one-step replica symmetry breaking (RSB) is used within the static ansatz. 
This formalism has already been applied in a theory in which m (Parisi's 
block-size parameter) is taken as a constant 13 . Now, it is extended to 
consider m as a variational parameter. In this case, the results show that 
RSB is present when T — > 0, in which the system is driven by quantum 
fluctuations. 

The infinite range Ising spin glass (ISG) in a transverse field has been ex- 
tensively studied as a theoretical quantum spin glass model. The strength of 
quantum fluctuations is adjusted by the transverse magnetic field. An experi- 
mental realization of this model is the quantum spin glass system LiHo x Yi^ x F^ 
[J. Moreover, there are now several disordered Cerium and Uranium alloys |2j 
in which it is possible to find a spin glass phase (SG) as well as a Kondo behav- 
ior around a quantum critical point (QCP). It is not our intention to discuss 
this particular complicated problem. However, a proper fermionic formulation 
of the SG problem in a transverse field can also be useful for this mentioned 
class of strongly correlated problems [3]. 

Several techniques have been used to treat the ISG in a transverse field. 
Nevertheless, they show some controversial results concerning the local stability 
of replica symmetry solution jUEHHIIIllBlE]- For instance, it is an open question 
whether or not the quantum tunneling through the barriers between the many 
degenerated thermodynamic states in the free energy landscape is able to restore 
the replica symmetry (RS). The quantum ISG model treated by the Trotter- 
Suzuki formalism in the static approximation 10 has shown stable RS solutions 
in a small region close to freezing temperature T c pi]. However, the replica 
symmetry breaking (RSB) is found in whole SG phase when the Trotter-Suzuki 
formalism is used without the static approximation, in which the dynamic spin 
self-interaction is treated by a numerical calculation [S]. For this case, the 
ordered SG phase has been investigated by one-step RSB solution, (i 7 in which 
it is suggested that the quantum tunneling does not restore the RS. On the other 
hand, this quantum problem has also been investigated by using imaginary-time 
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replica formalism under the static ansatz. This approach has shown that the 
RS quantum SG phase is stable in almost whole SG phase [Sj. Furthermore, the 
SG quantum rotors as well as the transverse field ISG have been studied within 
a Landau theory in which the dynamic effects are considered in an analytical 
approximation near the critical temperature In this treatment, the RSB 

occurs in the SG phase at finite temperatures and it suggests suppression of the 
RSB when T — > 0. 

Recently, the fermionic ISG in a transverse field Y has been studied using 
functional integral formalism with Grassmann variables 9 . This problem has 
been treated within the static approximation. The phase diagram found in Ref. 
9 shows that the T c decreases toward a QCP when T reaches a critical value. In 
this formalism, the RS solution is unstable in whole SG phase 9 in disagreement 
with Ref. [HI- In a subsequent work still using the static approximation, the 
replica symmetry has been broken within the one-step RSB like-Parisi's scheme 
[T2| . Nonetheless, the block -size parameter m has been taken as a temperature 
independent ^H] ■ The results within this particular RSB treatment suggest that 
the RS is restored at T = when Y > However, it arises some questions 
whether m is taken as a variational parameter as in the original Parisi's scheme 
[T^] for r > 0. Are the results using the original Parisi's scheme of RSB changed 
significantly when compared with those obtained in Ref. for T > 0? Are 
the quantum effects able to restore the RS in this fermionic formulation? 

The purpose of the present work is to investigate the one-step RSB in the 
fermionic ISG in a transverse field using the original Parisi's scheme. In particu- 
lar, we compare it with the alternative procedure of the same problem proposed 
in Ref. [13] ■ The problem is formulated in a path integral formalism as in Ref. 
1 1 H j . The static approximation is used to treat the spin-spin correlation func- 
tions. However, the parameter m is taken as a saddle-point equation, which is 
the main difference between this work and Ref. J3| . In our quantum spin glass 
treatment, the static approximation is an important point. It is recognized that 
the static approximation can yield inaccurate quantitative results at very low 
temperature. However, because of a lack of a suitable analytical method to go 
beyond the static ansatz in whole SG phase in this fermionic formalism, there 
are still several quite recent papers studying the ISG in transverse field within 
this particular level of approximation [5] This paper reports results for 

the one-step RSB order parameters, free energy and entropy. Our results agree 
partially with those of Ref. ^Sj- Nevertheless, the RS is not restored at T — > 0, 
which is in agreement with Refs. [HIE]. 

The infinite range ISG in a transverse magnetic field Y is described by the 
Hamiltonian: 

H=-Y / J ij S?S*-2Yj2S? (1) 

where the sums are run over the N sites of a lattice. The exchange interaction JV,- 
among all pairs of spins is a random variable that follows a Gaussian probability 
distribution with mean equal to zero and variance 16 J 2 . The spin operators in 
Eq. ||TJ are defined in terms of fermion operators [5]: Sf = \{n^ — rijjj, Sf — 
|[ c iT Ci i + c Il Ci T ] where h% a — c\ a Ci a gives the number of fermions at site i with 
spin projection a =t or J.. c| and Ci a are respectively the fermion creation and 
annihilation operators. 

The Hamiltonian Q is defined on a fermionic space where the operator S z 
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has two magnetic eigenvalues and two non-magnetic eigenvalues (empty and 
double occupied site). Therefore, there are two states that do not belong to the 
usual spin space. However, it is possible to express the partition function of a 
spin system (2S-model) in terms of the partition function of a corresponding 
fermionic system by using a restriction that considers only sites occupied by 
one fermion |15|. This work considers two models: the 4S-model without any 
restriction about the non-magnetic states and the 2S-model in which the non- 
magnetic states are forbidden. The restriction in the 2S-model is obtained by 
using the Kronecker delta function |15| . Therefore, the partition functions for 
both models are written in a compact form as 0: 

Z{y} = J D{4>*4>) fl i- JJ d Xj e-^e A ^ (2) 

with the action 

A{y} = jf dr{Y^ + Vj]4>i*(r) 



3,a 

-F(^ CT (r),^(r))} 



(3) 



/3 = 1/T (T is the temperature), yj = ixj for the 25- model or yj — for the 
45- model, which corresponds to the half-filling situation. 

In Eq. the Fourier decomposition of the time-dependent quantities is 
employed. The configurational averaged free energy per site is obtained by using 
the replica method: (3F = —-k lim ((Z{y} n ) j, . — l)/n where the replicated 

partition function is 

J 1 a, 7 ( 4 ) 

+ lnA{y})] 



with J dU = Y\ a dq a y(/3Jy / N/2iT) where a and 7 are replica indices 
running from 1 to n. It is assumed the static approximation, therefore: 



A{2/} = ny o -2^ e ~ Va J ^[C^]exp[i?], (5) 

s=j2^dr+^ 2 J 2 Erf^ ( 6 ) 

a 0,7 

with the definitions: 

5 * = 5 ( g ) 

^(w) = f^laC^) 0*a(' 4 Oj * s a Grassmann spinor, cr" and <j_ z are the Pauli 
matrices, and u denotes fermionic Matsubara frequencies. 
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The set of integrals in J dll has been exactly performed in the thermody- 
namic limit by the steepest descent method. The simplest conjecture to the 
replica matrix {Q} for overlaps in the spin space is the replica symmetry ansatz 
(RS). However, it produces an unstable solution in the ordered spin glass phase 
An attempt to solve this problem is the so called method of replica sym- 
metry breaking (RSB) |12| . 

In the Parisi's scheme of one-step RSB (1S-RSB), the replica matrix {Q} is 
parameterized dividing the n replicas into n/m groups with m replicas in each 
one, such as: q aa = r and 

q ay = qiifl(a/m)=l(j/m) 
q ai = qa if I{a/m) ^ I(j/m) 

where I(x) gives the smallest integer which is greater than or equal to x. The 
parameter r is the replica-diagonal spin-spin correlation r = (S a S a ). The 
physical meaning of m is that two different spin glass order parameters qo and 
qi are mixed up in the ratio m : 1 — to. 

The parametrization JHJ) is used to sum over the replica index. It produces 
quadratic terms in Eq. JHJ that are linearized introducing new auxiliary fields 
in Eq. ©. The functional Grassmann integral is now performed and the sum 
over the Matsubara frequencies can be evaluated like references HO EH- Finally, 
the restriction condition over the number of states admitted in each model is 
explicitly adopted. This procedure results in the following expression to the free 
energy: 

[3F s = Pf[m{ql-ql)+r*-ql] 

1 f> oo poo 

(10) 

/ Dzhx l Dv{2K s ) m 

m J -OO J -QO 

with 

K s = —— + / £>£cosh(/?</A) (11) 



where s (= 2 or 4) represents the number of states admitted in each model, Dx = 
dx e- x2 / 2 /V^ (x = z,v or £), A = h 2 + T 2 , h = J{^/2q^z + ^/2{ qi - q )v + 
y/2(r — <?i)£) and the parameters qo, qi, r and m are given by the extreme 
condition of the free energy ifTrljl . 

The set of saddle-point equations for the parameters q$, q\, r and m can 
numerically be solved. Here, the equation for m is self-consistently solved with 
the parameters qo, q\ and r |12| instead of Ref. |13j where m is taken as a 
constant value belong to the interval < m < 1. 

The 2S model recovers the classical results for 1S-RSB 12J when T = 0. 
When r > 0, the replica-diagonal correlation r becomes temperature dependent 
and r must be evaluated thogether with the other parameters (qa, q\, to). 
The results are shown in Fig. ^ The parameters qo, qi, m and S = q\ — qo 
(which gives 1S-RSB parameter) are zero at T > T c where the replica symmetry 
solution is valid. The transverse field reduces T c and the order parameters as 
well. The thermodynamic quantities are less sensitive on variations of m at T 
near T c than at lower temperatures. For example, the results obtained using 
the original Parisi's scheme show that 6 increases when T decreases from T c 
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: f = l 




Figure 1: Parameters qo, qi, 5 = qi — qo, r and m as a function of T/J for 
the 2S model. The upper, the middle and the lower panels are, respectively, for 
T/J = 0.0, 0.5 and 1.0. The dotted lines show the RS solution for r. 

in agreement with Ref. |13| . This occurs because right below the transition 
temperature the remaining RSB order parameters have very weak dependence 
on the solution of parameter m. On the other hand, at lower temperature, 
the original Parasi's scheme shows that the parameter m becomes increasingly 
important. For instance, the parameter S is greater than zero and the difference 
between r and r^s (for T > 0) is increased when T decreases by contrast with 
Ref. where S becomes zero when T — > 0. Therefore, our results indicate 
that the validity range in temperature of RSB procedure proposed in Ref. ^Sj is 
restrained to a region quite close to the critical temperature T c . As consequence, 
any conclusion about whether or not the RS is restored in Ref. ^3] is limited 
by the procedure itself. On the contrary, our results obtained by the original 
Parisi's scheme (in the same static approximation) suggest that the RSB effects 
are present even at T — > as it can be seen by the extrapolation of curves in 
Fig. 1 to T -> 0. 
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Figure 2: Parameters q$, q\ 1 5 = q\ — qo, r and m as a function of T/J for the 
4S model. It is used the same convention as Fig. 1 for the dotted lines and 
values of T. 

Fig. Elshows the order parameters for the 4S model. The behavior of qo, q±, 
6 and m is qualitatively the same as that one of the 2S model. However, r is 
temperature dependent for all values of T due to the presence of non-magnetic 
states. On the other hand, the ground state energy of the 4S model when T = 
is the same as that of the classical 2S model .16.. The results of the 2S and 
the 4S models are converging to the same values when T — > even for T > 0. 
Therefore, it could be expected that in the ground state of the 4S model every 
site carry a single fermion. Hence the empty and the double occupied sites do 
not contribute to the thermodynamic quantities at T — > 0. 

The free energies for 2S and 4S models when T = 0.5 J and 1.0 J are exhibited 
in Fig. |21 The 1S-RSB can be compared with the RS solution. The difference 
between these two approaches increases when T decreases from T c . The 1S-RSB 
ansatz gives higher values for the free energy, and its entropy (S — —dF/dT) 
is positive for the range of temperatures analyzed in Fig. [21 The increase of T 
reduces the difference between the RS and the 1S-RSB. It could be conjectured 
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Figure 3: Free energy as a function of T/J for the 2S and 4S models. The upper 
and the lower panels are, respectively, for T/J — 0.5 and 1.0. The solid lines 
are results for the 1S-RSB, while the dotted lines are for the RS approach. The 
insert shows the entropy at lower temperatures. 



that when the strength of Y is enhanced quantum fluctuations assume a relevant 
role inducing tunneling between the metastable spin glass valleys. Thereby, it 
could reduce the importance of the replica symmetry breaking. Nevertheless, it 
is still not able to restore the RS solution. Although these results are limited by 
the static ansatz that is expected to yield inaccurate quantitative results at very 
low temperature, they are in good qualitative agreement with those obtained by 
the Trotter-Suzuki formalism without the static ansatz O E| when T decreases 
from T c and by the Landau theory of quantum spin glass rotors ^T] when T 
increases towards the QCP. 

To conclude, it is studied the fermionic Ising spin glass model in a transverse 
field r within a Grassmann path integral formalism. The static ansatz and the 
replica trick with one-step replica symmetry breaking scheme are used. Recently, 
this problem has been analyzed by the same formalism and the same static 
approximations, but within a theory in which the Parisi block size parameter < 
m < 1 is taken as an independent parameter of temperature . This method 
in Ref. for finite and fixed values of m leads to the conclusion that the RSB 
parameter decreases towards zero at low temperature and the RS is restored at 
T = when T > 0. We have compared the procedure proposed in Ref. ^3] 
with the original Parisi's scheme which indicates that this alternative procedure 
is reliable mainly in the vicinity of the critical temperature. Therefore, any 
extrapolation to lower temperature is limited by the procedure itself besides the 
static approximation. Moreover, when quantum fluctuations become important, 
our results provide evidences that the original Parisi's RSB scheme should be 
used to improve the physical description of the ordered spin glass phase treated 
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by the replica method. 
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